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ABSTRACT 


Several flight-critical aerodynamic problems such as buffet, flutter, stall, and wing rock are strongly 
affected or caused by abrupt changes in unsteady aerodynamic loads and moments. Advanced sensing 
and flow diagnostic techniques have made possible simultaneous identification and tracking, in real-time, 
of the critical surface, viscosity-related aerodynamic phenomena under both steady and unsteady flight 
conditions. The wind tunnel study reported here correlates surface hot-film measurements of leading 
edge stagnation point and separation point, with unsteady aerodynamic loads on a NACA 0015 airfoil. 
Lift predicted from the correlation model matches lift obtained from pressure sensors for an airfoil 
undergoing harmonic pitchup and pitchdown motions. An analytical model was developed that 
demonstrates expected stall trends for pitchup and pitchdown motions. This report demonstrates an 
ability to obtain unsteady aerodynamic loads in real-time, which could lead to advances in air vehicle 
safety, performance, ride-quality, control, and health management. 
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x/c location of SEP 

angle-of- attack, degree 

time derivative of angle of attack 

relaxation time constant 

total time delay of SEP and reattachment 



INTRODUCTION 


In the past several decades, the flight envelope has been expanding for military aircraft into high 
angles of attack and increased angular rates. Entry into increased angle-of-attack (a) regimes has 
consequently shifted the emphasis from steady-state, to quasisteady , 1 to unsteady aerodynamics. With 
this shift, numerous complex problems associated with unsteady aerodynamics have come to the 
forefront. Problems such as flutter, buffet, abrupt wing stall, and wing rock adversely affect 
maneuverability and controllability of the aircraft. Although recent interest in solving unsteady 
aerodynamic problems has been spurred by military aircraft, civilian aircraft have had similar problems 
when subjected to adverse flight conditions such as crosswinds, gusts, and vortex encounters. 

The shift from quasisteady to unsteady aerodynamics lies partly in the inclusion of viscous effects. 
Viscous effects such as boundary layer transition from laminar to turbulent flow, flow separation and 
reattachment, and vortex dynamics play a large role in the performance, safety, and control of aerospace 
vehicles. For each of the unsteady aerodynamic problems mentioned earlier, there are associated viscous 
phenomena causing it: buffet can be caused by separated flow, inducing irregular oscillations to aircraft 
components; abrupt wing stall is caused by abrupt asymmetric flow separation on the wing, causing loss 
of lift on one wing; flutter can be caused by flow separation or shock impingement; and wing rock can be 
caused by asymmetric leading-edge vortices, inducing a limit-cycle oscillation in roll. 

Traditionally, viscous effects have been very difficult to measure, yet arguably are the most untapped 
source of information for aerodynamics study. Large-scale flow structures emanate from the boundary 
layer. Observing the off-surface flow phenomena is useful in a qualitative manner, but probing the 
phenomena as it develops in the boundary layer provides the opportunity to understand fundamental flow 
physics as well as the ability to profoundly influence the resulting flow structures through small-scale 
modifications to the boundary layer. 

Conventionally, instruments such as strain gages, accelerometers, pressure transducers, pressure 
sensitive paint, and even microphones are used to measure flow state and relevant aerodynamics. Each 
has distinct disadvantages. Strain gages and accelerometers indirectly derive the aerodynamic forcing 
function by measuring structural response. For load balance measurements, structural response of the 
system is only as accurate as the model used to map the measured strains to the imposed loads. 
Pressure-based instrumentation is typically challenging to set up and is subject to pneumatic latencies. 
Although pressure sensitive paint provides continuous spanwise pressure information, it suffers the same 
problems as other pressure-based instrumentation: significant setup time and strict temperature 
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conditions. For unsteady wind tunnel tests, each of the above instrumentation types has major issues that 
limit its use. Temperature, pressure, and density are just a few of the environmental variables that change 
often in flight, so it would be beneficial to use a measurement technique that requires the minimal amount 
of calibration possible. 

Hot-film anemometry can measure viscous effects like shear stress. This same method can also be 
used to characterize viscous phenomena such as transition, flow separation and reattachment, turbulence, 
and shock impingement. Hot-film sensors are also nonintrusive and can be installed on bodies without 
major structural changes. However, conventional anemometry has certain limitations. The large time 
constant of hot films (as opposed to hot wires) makes constant current anemometry unsuitable for use 
with hot films . 4 Constant temperature anemometry can handle the large time constant but is complex to 
set up. An improvement upon conventional anemometry is the constant voltage anemometer (CVA). This 
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instrument has a simple setup and high sensitivity with large bandwidth, therefore, it is an excellent 
choice for this application. 

Another advantage of using hot-film anemometry for unsteady loads measurement is the ability to use 
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a flow diagnostic technique called phase-reversal signatures (PRS), that does not depend on the 
absolute values of shear stress, thereby significantly reducing the need for calibration. Monitoring the 
PRS with the CVA enables quantitative in-flight measurement of the flow state in real-time. 

This report describes an approach that allows quantitative measurement and tracking of surface flow 
phenomena in order to obtain loads with minimal calibration requirements. This approach can even be 
used in the presence of changing environmental variables such as temperature and pressure, since the 
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basic components of the technology behind the approach have been demonstrated in flight. This 
report presents a description of the wind tunnel setup (model, instrumentation, and data acquisition 
details), data analysis results, and interpretation of the unsteady aerodynamic data. 

THEORY 

The objective of this report is to make use of the relationship between surface flow phenomena, 
leading-edge stagnation point (LESP) and flow separation point (SEP), and aerodynamic loads. The 
important step then, is to develop that relationship empirically by correlating LESP and SEP with 
aerodynamic loads. The theory of obtaining a polynomial relationship as a correlating function of LESP 
and SEP to aerodynamic loads can be traced to a mathematical state-space representation by Goman and 
Khrabrov. 14 Goman and Khrabrov describe the mathematical state-space representation for the unsteady 
aerodynamic problem as: 

1. State Equation 

TjX + x = x 0 (a-x 2 a) (1) 

2. Output Equation 



where, Xq is the nonlinear, steady-state dependency of SEP on a; x is the instantaneous location of the 
SEP, T| is the relaxation time constant which defines the transient aerodynamic effects; x 2 is the total 
time delay of the SEP and reattachment caused by quasisteady aerodynamic effects, such as the 
circulation and boundary-layer convection lags. The driving function is the steady-state location of the 
SEP, which is generally a nonlinear function of a and usually estimated from wind tunnel tests. From the 
output equation, the lift coefficient is shown as a function of a and instantaneous SEP. 

In the above state-space model, for equation 1, there are two unknown time constants that need to be 
determined using parameter estimation techniques with a large dataset consisting of experimental 
unsteady aerodynamic data at various unsteady flow conditions. For each condition, equation 1 needs to 
be integrated to obtain a value for x. At the same condition, a and x are input to equation 2 which yields 
the lift coefficient. 
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This state-space method for determining unsteady loads has been validated by Goman and Khrabrov 
and Fan and Lutze 15 for various unsteady flow conditions involving flow separation and vortex 
breakdown using wind tunnel measurements on an airfoil, a flat delta wing, and an aircraft configuration. 
So, the state-space method allows the lift coefficient to be a function of only a and SEP. 

Inspired by Goman and Khrabrov’ s work, if a and SEP are measured in real-time, then the lift 
coefficient can be computed directly from equation 2. The ability to directly measure a and SEP would, 
therefore, eliminate the need for extensive and expensive wind tunnel testing required to obtain the time 
constants. If LESP can be substituted for geometric a, then hot-film sensors could be used to provide the 
necessary components for a polynomial approach. Hot-film sensors and the CVA can be used to provide 
real-time measurements of LESP and SEP. The LESP location is linearly proportional to geometric a in 
the prestall region. Near and after stall, LESP movement deviates from geometric a. Instead of using the 
output equation, a first-order polynomial approximation of lift as a function of LESP and SEP is used. 
This approach can be intuitively justified by observing the relationship of the leading-edge and rear 
stagnation points with circulation. According to the Kutta-Joukowski condition, 16 the circulation around 
an airfoil is determined by imposing the rear stagnation point to be located at the sharp trailing edge. 
When the rear stagnation moves upstream on the suction side because of flow separation, there is a 
significant loss in circulation which results in decreased lift and increased profile drag. So, intuitively, the 
LESP and SEP should be well correlated with circulation. 

For this approach to work, the steady-state dependency of LESP and SEP on a needs to be measured 
to create a polynomial that correlates steady-state loads with the LESP and SEP flow features. Once the 
first-order polynomial is obtained, LESP and SEP data can be input into the polynomial equation to 
extrapolate the unsteady lift coefficient (C L ). If the unsteady lift from the polynomial matches that from 
the pressure sensors then the method to measure loads from LESP and SEP is validated. 

TECHNICAL APPROACH 

Multiple hot- film sensors and the CVA were used to identify the LESP and SEP on a fixed and 
oscillatory NACA 0015 airfoil. A three-step approach was used in correlating LESP and SEP to unsteady 
loads on a NACA 0015 airfoil. The first step was to capture the flow features quantitatively through the 
airfoil stall region at steady- state and at different positive and negative pitch rates. The second step was to 
correlate the LESP and SEP with measured lift. The third and final step was to obtain a polynomial 
relationship between the LESP and SEP and the measured lift such that all observed unsteady cases are 
functions of the viscous phenomena. The United States Air Force Academy (USAFA) subsonic wind 
tunnel was utilized to set up the instrumentation, conduct several unsteady experiments, and collect data 
related to the viscous phenomena and lift. 

Model and Instrumentation Description 

The model tested was a symmetric 6-inch chord NACA 0015 airfoil. The airfoil was instrumented 
with 14 pressure ports and 104 hot-film sensors. Along one side of the airfoil, there were 14 pressure 
ports from approximately 2-percent to 90-percent chord and 70 hot-film sensors from approximately 
1 -percent to 97-percent chord. On the other side, the airfoil was instrumented with hot-film sensors 71 
through 104 from approximately 1 -percent to 40-percent chord. The leading edge is between sensors 70 
and 71. Figure 1 shows the instrumented airfoil within the test section. 
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Figure 1 . Instrumented airfoil in United States Air Force Academy subsonic wind tunnel. 


An oscillating mechanism was constructed for the test. A proportional integral derivative (PID) 
controller was made to control the airfoil motion for accurate pitch rates. The angles were tracked using a 
quadrature encoder. 

Once the instrumented NACA 0015 airfoil was integrated with the oscillating mechanism in the wind 
tunnel, the remaining associated instrumentation was connected. All 104 channels of hot-film sensors 
were each connected to a dedicated CVA. The CVA shear stress output was recorded onto analog tapes, 
along with the 14 pressure port data and the quadrature-encoded pitch angle. The hot- film, pressure, and 
angular data were recorded simultaneously onto analog tape and later digitized at a sampling frequency 
of 6 kHz and 12 bit analog-to-digital conversion. 

Experimental Approach 

The first step in the experimental approach was to capture the LESP and SEP quantitatively through 
the airfoil stall region at steady-state and at different positive and negative pitch rates. 

Wind tunnel testing was performed at atmospheric conditions from speeds of 50 ft/sec to 100 ft/sec. 
Data were obtained over the a range of minus 40 deg a to plus 40 deg a. For the steady-state cases, the 
entire range of angles was stepped through at one-degree increments. Since there are only 14 pressure 
sensors on one side of the airfoil, the symmetry of the airfoil had to be utilized for determining the 
pressure profile for the upper and lower surfaces. For the positive angles of attack, the pressure sensors 
provided the pressure profile for the upper surface, and for the negative angles of attack the same pressure 
sensors provided the pressure profile for the lower surface. 
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For the unsteady cases, the oscillation mechanism was sinusoidally driven: 


a(t) = A^l - cos(2jt/ t)\ (3) 

where geometric a is a function of time t, A is the amplitude, and /is the frequency. Three different 
frequencies were chosen and several amplitudes were used. The information in this report is limited to 
two frequencies, 0.4 Hz and 0.8 Hz, and one amplitude, 40 deg a. One cycle consists of the oscillation 
mechanism starting at 0 deg, rising sinusoidally to 40 deg, and then dropping back to 0 deg a. For the 
pressure sensors, the cycle must be mirrored for the other side as well, so that a goes from 0 deg to minus 
40 deg back to 0 deg. 

The LESP and SEP were identified primarily using PRS. However, as a redundant technique, the 
viscous phenomena were also identified by regions of minimum shear stress and the presence of a large 
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root mean square (RMS) signal nearby. 

After the LESP and SEP had been identified and tracked, they were correlated to the measured lift. 
Lift was calculated by integrating the surface pressures on one side of the airfoil twice: once while at a 
positive a and then again at the same a, except negative. The same technique was used to calculate lift 
from the unsteady pressure measurements. No latency corrections were made for the pressure 
measurements. 


RESULTS 

First, a typical shear stress plot is shown to explain how the different flow features are represented in 
shear stress measurement. Then, graphs tracing the locations of LESP and SEP are shown, with combined 
plots of lift measured from integrating pressures along the airfoil. These graphs show how the LESP and 
SEP are correlated with measured lift. Using the LESP and SEP as input to a model, lift is calculated and 
compared with measured lift. Also, shear stress plots are shown for different pitch rates to exhibit 
dynamic stall. 


Labeled Shear Stress Plot 

Figure 2 illustrates an example pointing out the different flow features to be tracked and analyzed. 
The shear stress output in figure 2 is from the airfoil, as described earlier, instrumented with hot-film 
sensors on both the suction and pressure side, undergoing a pitchup then pitchdown motion from 0 deg to 
40 deg and back to 0 deg a. High shear stress regions are colored red and low shear stress regions are 
colored blue. Flow bifurcation points such as LESP and SEP exhibit low shear stress. 
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Figure 2. Shear stress distribution labeled with significant flow features (100 ft/sec, 0.4 Hz). 


Focusing on the pressure side of the airfoil, observe the dark blue line moving initially from 0-percent 
chord at 0 deg to approximately 10-percent chord at 40 deg and then back to 0-percent chord at 0 deg. 
The differentiated line represents the instantaneous location and movement of the LESP. Now focusing 
on the suction side of the airfoil, observe another blue, very low shear stress line, moving initially from 
40-percent chord at 4 to 5 deg to 0-percent chord at 40 deg and then back to 40-percent chord at 3 to 
4 deg. This line is the SEP movement. As labeled on the plot, the patches of yellow on the suction side 
represent laminar separated flow. 

In figure 2 there are several trends to note concerning LESP and separated flow. For the positive pitch 
rate, the LESP moves away from the leading-edge linearly proportional to the increase in a to about 
10 deg and then begins to recede around 15 to 20 deg. Stall occurs at approximately 18 deg angle of 
attack. For the negative pitch rate, the LESP recedes at around 10 deg angle of attack. When stall occurs, 
the LESP is close to receding or moving back towards the leading-edge. Compare when stall occurs for 
positive and negative pitch rates. Notice how stall occurs at higher angles of attack for the positive pitch 
rate than for the negative pitch rate. 

With regard to the laminar separated flow, observe the relation between stall and the amount of 
laminar separated flow on the suction side. As a increases beyond stall, the SEP, darkest blue region on 
the suction side, moves toward the 0-percent chord resulting in leading-edge separation. The attached 
reversed flow from flow separation decreases rapidly beyond stall. Comparing the laminar separated flow 
for positive and negative pitch rate, at a positive pitch rate, the onset of leading-edge separation is 
delayed, and at a negative pitch rate, flow attachment is delayed. 


Stagnation and Separation Point Movement 

Although the LESP and SEP can be readily identified in a qualitative manner in the shear stress plots, 
for the purposes of this report it is necessary to trace their motions precisely. This tracing is done using 
phase-reversal signatures (PRS). In figure 3(a), we can observe a 180-deg phase-shift between sensors 67 
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and 68 and the same for 89 and 91. The LESP is between 3.3- and 4.1-percent chord (sensors 67 and 68). 
The SEP is very close to sensor 90 which is at 18.2-percent chord, between sensors 89 and 91 (16.7- and 
19.8-percent chord). Although it is not obvious from figure 3(a), there is a 180-deg phase-shift between 
sensors 90 and 91. This can be clarified by calculating the cross-correlation coefficient for sensors 90 and 
91 in figure 3(b). In figure 3(b), at zero lag, the correlation coefficient is negative, which is indicative of a 
180-deg phase- shift between sensors 90 and 91. The figures trace LESP and SEP for 100 ft/sec and 
0.4 Hz. 
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(a) 180-deg phase-shift for stagnation and 
separation point. 


(b) Cross-correlation of adjacent sensors at 
separation point. 


Figure 3. Phase-reversal signatures (PRS): 180-deg phase-shift and cross-correlation. 


In figures 4(a) and 4(b), the LESP and SEP movements are plotted as chord locations as a function of 
the angle of attack. Lines are used for figure 4(a) to help show the trend. As compared to the steady-state 
LESP movement, the pitchup LESP movement overshoots and the pitchdown LESP movement 
undershoots. The steady-state LESP drops at close to 15-deg a. The pitchup LESP drops much later at 
near 18-deg a. The pitchdown LESP drops near 9-deg a. Since the LESP is related to circulatory lift, it 
exhibits properties associated with dynamic stall. 
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(a) LESP movement. 


(b) SEP movement. 


Figure 4. LESP and SEP movement for steady-state, pitchup, and pitchdown (100 ft/sec, 0.4 Hz). 


The SEP movement shown in figure 4(b) shows how long flow takes to reattach for pitchdown. The 
pitchup SEP takes much longer to get to the leading edge, longer than even the steady-state SEP. So, 
leading-edge separation is delayed for the airfoil pitching up. 

Figures 5(a) and 5(b) look very similar. Figure 5(a) is lift calculated from measured pressures for 
pitchup and pitchdown motions. Ignoring the pneumatic lag in the pressure output and the measured lift 
coefficient, one can see, by comparing figures 5(a) and 5(b), the clear correlation between the hysteresis 
in lift and the corresponding hysteresis in LESP movement. 



0 5 10 15 20 25 30 


Angle of attack, a, deg 


040157 


• Pitchup 
T Pitchdown 










tI*' 

r* 



#* 


T 




• 

• • 


T 

f 



▼I 1 1 

,• 





T* 1 

▼ 

F 





¥• 







0 5 10 15 20 25 30 


Angle of attack, a, deg 


040158 


(a) Lift measured from integrating pressures. 


(b) LESP movement. 


Figure 5. Lift and LESP movement for steady-state, pitchup, and pitchdown (100 ft/sec, 0.4 Hz). 
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Plotting the measured lift on top of the LESP movement clearly shows the correlation. Figure 6 
shows the lift and LESP movement are plotted together. Figure 6(a) shows that the steady-state LESP 
movement follows lift linearly to the stall angle. Before stall, for an increase in lift coefficient of 0.2, the 
LESP moves approximately 1 -percent chord away from the leading edge. The pitchup and pitchdown 
LESP movements in figures 6(b) and 6(c) have similar trends to the pressure-integrated lift. The next step 
is to use the LESP and SEP as inputs into a polynomial function that approximates the lift coefficient. 
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(a) Steady-state. 



(b) Positive pitch rate (0.4 Hz). 


Figure 6. LESP movement and lift (100 ft/sec). 
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(c) Negative pitch rate (0.4 Hz). 


Figure 6. Concluded. 


Lift Polynomial 

The following equation is a first-order polynomial of the lift coefficient as a function of LESP and 
SEP x/c locations: 


C L( x lesp’ x sep) = a 0 x lesp + a \ x sep ( 4 ) 

where, the a coefficients need to be determined through a least-squares fit of steady-state input: LESP, 
SEP, and lift. A distinct prestall and poststall lift distribution was observed. Therefore, two polynomials 
were needed to describe the lift curve. Stall is empirically determined to be when SEP reaches 5-percent 
chord. Figure 7 shows both the resulting prestall and poststall polynomial equations for 100 ft/sec. The 
prestall equation is: 


C L{ X lesp’ X sep) = 0205lx lesp + °- 0002x sep 


( 5 ) 


and the poststall equation is: 


C L\ x lesp ,X sep) = 0AU3x lesp + 0M06x sep 


(6) 
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Figure 7. Comparison of pressure (solid lines) to polynomial fit (symbols) (100 ft/sec, 0.4 Hz). 


To predict lift, the locations of LESP and SEP in the unsteady pitchup and pitchdown cases are then 
inserted into the above equations, as xj esp and x sep respectively (fig. 7). The unsteady lift output in 
figure 7 clearly exhibits the same hysteresis as the lift curves. The predicted lift curves are not very far 
from the pressure-based lift coefficients. As a first cut, flow features such as LESP and SEP seem to 
represent the lift curves. Observe that for the prestall the polynomial equation, a j = 0.0002, is very 
small compared to a Q = 0.2051 . This observation shows that LESP has far more influence, in a 
least-squares sense, than SEP has on lift in the prestall angles of attack. Therefore, the LESP by itself can 
approximate the lift curves fairly well before stall occurs. The location at which LESP recedes towards 
the leading edge seems to be a precursor to stall. We shall use this observation in the next section when 
comparing increasing pitch rates. 


Increasing Pitch Rates 

The next issue is whether the flow features are consistent with published results as regards dynamic 
stall. At high positive pitch rates, increased lift is expected, and, at high negative pitch rates, decreased 
lift is expected. To verify whether or not the flow features are able to produce the expected trends for 
different pitch rates, two shear stress plots at two different frequencies are shown and analyzed. 

Figures 8(a) and 8(b) have labels on the left showing approximately where the LESP starts to recede 
for pitchup and pitchdown motions. For the positive pitch rates, the LESP recedes around 18 deg a for 
0.4 Hz and 19 deg a for 0.8 Hz. So, for positive pitch rates, increasing the pitch rate, increases the angle 
at which LESP recedes. For the negative pitch rates, the LESP recedes near 9 deg a for 0.4 Hz and 
8 deg a for 0.8 Hz. So, for negative pitch rates, increasing the pitch rate, decreases the angle at which 
LESP recedes. If we were to use the recession of the LESP as a precursor to stall, then at higher positive 
pitch rates, stall is delayed. Similarly, at high negative pitch rates, the airfoil stayed under stall conditions 
longer. For positive pitch rates, as the frequency increased, the flow was attached longer. For the negative 
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pitch rates, as the frequency increased, flow reattachment was delayed further. Therefore, the LESP 
movement provides the trends of dynamic stall. 



(a) Pitch rate (0.4 Hz). 


(b) Pitch rate (0.8 Hz). 


Figure 8. Shear stress profile: Increasing pitch rates (100 ft/sec, 0.4 Hz and 0.8 Hz). 


CONCLUDING REMARKS 


For a pitching NACA 0015 airfoil, leading-edge stagnation point (FESP) and flow separation point 
(SEP) movement was tracked using both shear stress plots and phase-reversal signatures (PRS). A 
steady-state polynomial relationship between lift and the locations of FESP and SEP was empirically 
developed. Using the steady-state polynomial relationship and the locations of LESP and SEP at airfoil 
pitchup and pitchdown, unsteady lift was extrapolated and compared to lift calculated from pressure 
sensors. This comparison showed the lift derived from the LESP and SEP tracked the lift calculated from 
pressure sensors. The LESP movement was further shown as similar to lift curves. The recession of LESP 
towards the leading edge was indicated to be a precursor to stall. The LESP recession at different pitch 
rates showed trends similar to that of dynamic stall. Such data indicate it is possible to obtain unsteady 
loads in real-time from LESP and SEP. These results have been obtained in incompressible flow. Further 
research needs to be conducted to consider compressible and Reynolds number effects. 

Identifying and tracking flow features such as LESP and SEP in real-time provides a clear benefit in 
the unsteady aerodynamics field. Analytically, the tracked LESP and SEP can provide input to an 
empirical code that more accurately predicts unsteady effects. For wind tunnel and flight test studies, 
tracking LESP and SEP can be a way to capture any scale effects since the LESP and SEP can be 
accurately tracked in the wind tunnel and in-flight without any calibration. Any change in unsteady 
viscous behavior resulting from scale differences between the model and actual flight vehicle exhibits 
itself through LESP and SEP. Such changes can be recorded, analyzed, and incorporated back into 
empirical codes. 
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For flight experiments, tracking the LESP and SEP can allow for surface flow diagnostics that help 
understand and prevent abrupt wing stall, wing rock, onset of flutter, and buffet. Incorporating real-time 
flow diagnostics as an integrated smart sensor would enable the development of a closed-loop active flow 
control solution, allowing for new levels of performance and safety for both military and civilian aircraft. 


Dryden Flight Research Center 
National Aeronautics and Space Administration 
Edwards, California 
June 1, 2004 
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